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Abstract 

It is shown that, in QCD, the same universal function r cnsp (i?, a s ) determines the in- 
frared behaviour of the on-shell quark form factor, the velocity-dependent anomalous dimen- 
sion in the heavy quark effective field theory (HQET) and the renormalization properties of 
the vacuum averaged Wilson lines with a cusp. It is demonstrated that a combined use of 
the methods developed in the relevant different branches of quantum field theory essentially 
facilitates the all-order study of the asymptotic and analytic properties of this function. 
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1. Introduction 



A considerable effort is being made now to construct a new approach to the heavy-quark physics. The basic 
idea [1, 2] is to use the fact that the masses of the heavy quarks c and b are essentially larger than the QCD 
scale A, and to start (when appropriate) with an effective field theory [1, 2], in which the heavy quark masses 
are infinite: tuq — > oo. 

In the infinite mass limit, the components of the heavy quark momentum p^ tend to infinity together with the 
quark mass: p <~ toq, the ratio p^/mQ approaching a fixed value v^. The 4- velocity w M = (1/vl — v 2 , v/ yl — v 2 ) 
just specifies the direction in which the infinitely heavy quark is moving, and the heavy quark effective theory 
(HQET) [3, 4] has a peculiar property that each effective quark field h v (x) is characterized by a 4- velocity vector 
V 

The electroweak interactions can change the velocity of the heavy quark v 1 — > v 2 and, as noted by Isgur and 
Wise [2] , it is the 4- velocity change rather than the momentum transfer that is an important dynamic variable 
for the heavy meson form factors, both for elastic (like j*D — > D and — ► B) and for weak (B — > De^-type) 
transition form factors. In other words, the QCD currents Q2^Qi are substituted in the HQET by the effective 
currents h V2 Th Vl containing fields related to the two velocities v\ and x>2- 

An important observation made in [5] is that the effective currents h V2 Th Vl have a universal (the same for 
all T- matrices) velocity-dependent anomalous dimension j w (a s ). The velocity change is characterized by the 
scalar product w = (viv 2 ) or, since the velocity vectors have a unit "length", v\ = v\ = 1, by the angle •& 
between them: w = cosh i9. 

As a matter of fact, the authors of [5] arc not the first who encountered an angle-dependent anomalous 
dimension in their calculation. Similar things happened before in two other and, at first sight, entirely different 
branches of quantum field theory. 

First, more than 10 years ago, Polyakov [6], trying to reformulate the dynamics of gauge fields in terms 
of string operators, has observed that the vacuum average of a Wilson loop having a cusp, possesses extra 
ultraviolet divergences. These specific divergences, "generated" by the cusp, are multiplicatively renormalizcd 
[7], and the relevant anomalous dimension T cusp ({), a s ) depends on the cusp angle d. 

The second example is related to the infrared (IR) behaviour of pcrturbative QCD which is known to be 
controlled by the subprocesses involving soft gluons. Since the soft gluons cannot change the momenta pi of the 
incoming and outgoing particles, one faces here a situation completely similar to that in the heavy-quark limit: 
the particles are just moving in fixed directions (specified by pi) through the gluonic cloud. The study of the 
IR behavior of perturbative QCD has not a purely academic interest: even if the IR singularities cancel, there 
remain finite IR-induccd contributions. In our papers [8, 9], we developed an approach that enables one to study 
these contributions using a renormalization group equation. The relevant anomalous dimensions Tm(&ij, a s ) 

depend just on the angles between the relevant momenta: cosh = (piPj)/ \J P 2 P 2 j ■ 

The IR-induced terms are especially important for high momentum transfer processes, when Q 2 ^ = ~{pi~ 
Pj) 2 p\ ~ p 2 = p 2 , i.e., when some angles tfy are large. We established that, in this limit, the anomalous 
dimension Tm^ij, a s ), to all orders in a s , is linear in or, what is the same, in log(Q 2 j /p 2 ): 

r Jfl (0y,a s ) = K{a s )\og{Ql/p 2 ) + A(a s ). (1) 

The function K(a s ) plays an important role in the IR-induced and the double logarithmic contributions 1 to 
different hard processes: the same combination 



K(a s ) = 




(where Cp =4/3 and Ca — 3 are quark and gluon color factors and Nf = 3 is the number of light quark 
flavours) governs the x — * 1 asymptotics of the deep inelastic structure functions [9] and the large perturbative 
corrections to the Drell-Yan cross section [11] and the e + e~— annihilation [12] . 

One might ask: are all these anomalous dimensions independent or there exists a connection between all 
of them? It is our goal in the present paper to study the interrelation between the infrared behaviour of 
perturbative QCD and the renormalization properties both of the Wilson loops (lines) and of the effective 

1 A function equivalent to K{a s ) is the basic object in another RG-type approach to the double- logarithmic 
effects developed by Collins [10] 
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heavy quark currents in the HQET, to demonstrate that the apparently unrelated functions, each playing a 
fundamental role in its field, are precisely identical and correspond to the QCD analogue of the bremsstrahlung 
function known from the early QED years. 



2. Infrared singularities and the heavy quark limit 

A simple but essential fact is that there are no velocity- or angle-dependent ultraviolet divergences in any 
standard QCD diagram. This means that the HQET has extra UV divergences: the heavy quark limit toq — > oo 
is singular. In particular, a logarithmic divergence in the HQET may result from a log(rng/A)-contribution in 
QCD, with A remaining finite in the heavy quark limit. However, since all the momenta p^ scale as mg in this 
limit: p^ — ► toqi^, the only combination that might be much smaller than Toq is (p 2 ~ ttiq), the heavy quark 
off-shellness, and the singularity may result only from a term like log(mg/(p 2 -nig)). Such a term is a standard 
infrared logarithm, with (p 2 — tUq) serving as an infrared cut-off. This indicates that the velocity-dependent 
divergences should be closely related to the IR behaviour of perturbative QCD. Indeed, studying the hadrons 
containing a single heavy quark, one is interested in the kinematic situation when the heavy quark energy differs 
from its mass by some finite amount E, which has a meaning of the average energy of the light quark(s), binding 
energy, etc. Thus, the virtuality of the heavy quark (p 2 — Toq) w IvhqE is much smaller than its mass squared: 

2 2 

A = «-2 « 1. 

m Q 

In the heavy quark limit, itiq — > oo, the parameter A tends to zero and, in this sense, the heavy quark is almost 
on-shell. 

It is well-known that, in gauge theories, the on-shell form factors have infrared (IR) singularities. The classic 
example is the elastic electron on-shell {p\ = p\ = to 2 ) form factor in QED. At the one- loop level one has, for 
the Dirac form factor 

a fm 2 \ ( Q 2 \ 
F(Q 2 , to 2 ) = 1 — — log ( -^2" j B ( — j j + regular terms) , 

where Q 2 = — (pi —pi) 2 is the transferred momentum, A is the (fictitious) photon mass serving as the IR cut-off 
parameter and B is the bremsstrahlung function 

^(S) = 4^io g | \ i-i. ( 3) 

Since A 2 is treated as a scale much smaller than to 2 , we have the same kinematic situation as in the heavy quark 
limit. 

The IR logarithms in QED are known to exponentiate when summed over all orders [13] 




F(()-.,n 2 ) = cxv{~ logl^B ( %^ ) } x (IB regular hrms). (4) 



In this expression, one can easily factorizc the "long distance" contribution (/i 2 /A 2 ) ^ B absorbing all the IR 
divergences, and the "short-distance" factor (m 2 /^ 2 ) ^ B . As a result, one can rewrite eq.(4) as 

where the IR regular factor C( m r, Qp) now has an explicit dependence on the IR factorization scale /i. 

The derivation of this factorized representation has a striking resemblance to the standard factorization 
procedure for the collinear mass singularities applied to deep inelastic scattering (in the form of the operator 
product expansion, OPE) and to other hard processes. Using it, one can factorize the original function, e.g., 
W n (Q 2 ,p 2 ), the n-th moment of the deep inelastic structure function, into the short-distance coefficient function 
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C n (Q 2 / jj, 2 ) and the long-distance sensitive matrix element fn(n 2 /p 2 ) absorbing all the collinear logarithms log p 2 
singular in the p 2 — > limit: 

W n (Q 2 , P 2 ) = f n (v 2 /p 2 )C n (Q 2 /v 2 ). 

The factorization scale fx separates the high- and low-momentum regions: the essential momenta k are 
k > \i for the C-function and k < \i for the /-function. An important fact is that the matrix elements of the 
local composite twist-2 operators O ni producing /„, possess extra UV divergences absent in the original theory 
(QCD). Of course, the factorization procedure automatically imposes the UV cut-off k < /i, and there are no 
actual infinities. But now one can make use of the fact that the ^-dependence of the /^-function is just the 
dependence on the UV cut-off parameter and incorporate the renormalization group to study it. That is how 
the n-dependent anomalous dimensions 7„(a s ) of the composite operators O n and/or z-dependent anomalous 
dimensions (evolution kernels) P(z,g) of the parton distribution functions f(x,[i 2 ) come into play. 

The analogy between the factorization properties of the QED formula (5) and the OPE-type factorization 
has a very deep reason [8]: in QCD (and in any similar gauge theory) one can factorize out the IR-singular 
terms of the on-shell quark form factor 

F(Q 2 ,m 2 ) — F IR x (IR regular terms). 

The IR-sensitive factor Fir accumulates all the effects due to interactions of the massive particle (heavy quark, 
if one studies the heavy quark limit) with soft quanta - gluons and massless quarks. The explicit expression for 
Fir is given by the vacuum average of the Wilson line operator W(C) 

Fir = (0|Tex P [ig £ dtv^A^v 2 t^j exp [ig J dt < A^vtt)) |0) = (0|TU(C)|0) (6) 

calculated along the classic path C of the massive ("heavy") particle. It goes from — oo along the direction of 
the initial momentum p\ = mvi to the interaction point and then along the direction of the final momentum 
P2 = mv2 to +oo. The soft gluons coupling directly to the heavy particle are described by the gauge potentials 
A,(x). ' 

To summarize, in this limit (one can treat it either as the IR limit A — > or as the heavy quark limit 
m — > oo), all the effects due to interactions of the massive particle with the gluons can be described by a Wilson 
line. This result has a natural interpretation: the wave function of a very massive particle propagating through 
the gluonic cloud acquires only a phase factor equal to the Wilson line. All other effects are suppressed by 
powers of A/m. In particular, they cannot produce contributions singular in the IR limit A — > 0. 

Just like in the OPE case, the factorization scale [i works as an ultraviolet cut-off for Fir, since only small 
momenta k < /j, contribute to it. And such a cut-off is really necessary because, as established by Polyakov [6], 
any Wilson loop with a cusp has extra UV divergences. 

3. Wilson lines and the heavy quark effective field theory 

The properties of the Wilson lines and loops where studied in detail at the beginning of the 80's. A very 
important result [14] is that one can interpret the vacuum average of a Wilson line W(C) as the propagator of 
one-dimensional fcrmions "living" on the integration path C: 

Tr W(C) = = J T>$(t)Thl>{t) Va(+oo)Va-oo)expHS) (7) 

where the action 5, in our case, is defined by 

/+oo q f° _ r+oo 

dti){t)—M)-g dt4>{t)v^A^(vit)ij{t)-g dt^v^A^v^it). (8) 
-oo at J-oo JO 

It is easy to show that the one-dimensional fermions defined in this way exactly coincide with the effective heavy 
quark field of the HQET: 

(0|/vM)|Mi) - (^(t)^(-oo)) = T exp (ig J dav^A^VKJ^j , t < (9) 
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and + 

(h,v 2 \h V2 (v 2 t)\0) = (ip(+oo)i>(t)) = Texp (ig j da v% A^a)^ , t>0. (10) 

The representation (7) is very convenient to analyze the general renormalization properties of the Wilson lines. 
In particular, it was found [7, 14] that for a smooth integration path all the UV divergences can be compensated 
by a renormalization of fields and vertices. 

The Wilson lines in the r.h.s. of eqs.(9) and (10) are not closed, and their end-points produce extra UV 
divergences. They can be multiplicatively renormalized, and the relevant anomalous dimension j en d(c(s) was 
calculated, in the Feynman gauge, to one loop in [7, 14] and to two loops in [15] 



From the relations (9) and (10) it follows that -)end{&s) coincides with the anomalous dimension of the effective 
quark field in the HQET and with that of the one-dimensional fermions in the theory defined by (8). The 
HQET calculation was performed in [16]. 

As mentioned earlier, each cusp of the integration path also generates extra UV divergences governed by the 
anomalous dimension T cusp that depends on the cusp angle i?. In the theory of one-dimensional fermions (8) one 
obtains T cusp as the anomalous dimension of the composite operator tp(+0)tjj(— 0). Using the relations (9) and 
(10) one can identify this operator with a heavy quark current h V2 (0)h Vl (0). Thus, the anomalous dimension of 
the heavy quark current in the HQET coincides with the cusp anomalous dimension of the relevant Wilson line 

T cusp (i),a s ) = j w (a s ), (11) 

provided, of course, that w = cosh <&. 

In general, it is a difficult problem to calculate Wilson loops, even in perturbation theory. However, it is 
possible to evaluate (0|W(C)|0) = Fir using the explicit form of the contour C defined by eq.(6). The path 
C has an infinite length and, as a consequence, (0|W(C)|0) has infrared divergences. Moreover, (0|IU(C)|0) 
depends on the directions of the rays forming C. Note now that the directions defining the line C are specified 
by the dimensionless ratios — p^/m and v% — pt^/m . Hence, the dimensionful parameters that explicitly 
appear in the vacuum average (6) are the IR cut-off A and the factorization scale /x, the latter, as emphasized 
earlier, working as an ultraviolet cut-off for the new divergences generated by the cusp of the path C at the 
point 0. This means that the UV scale (i and the IR scale A should basically appear in the ratio. There is also 
the ordinary renormalization scale p,R fixing the definition of the running coupling constant 5(/x_r): 



F IR = F IR ( f , — , (v lV2 ), g(p R )) = F IR (ft, — , , gfa) ) . 

A pr J \X pr m 2 



After taking p — pr, the IR factor Fjr, being equal to the vacuum average of a Wilson loop, obeys, in general, 
the RG equation 

where the anomalous dimension Tir controls the IR induced logarithms log(/z/A). Since the path C has a cusp, 
one identifies Tjr with the cusp anomalous dimension: 

Tir (^Us) = as). (12) 

Furthermore, since the IR singularities of Fir coincide with those of the original amplitude (e.g., with 
those present in the on-shell form factor), the cusp anomalous dimension T CUS p(-&, a s ) must reproduce the 
bremsstrahlung function (3): 

T cusp (#, a s ) = ^C F B (%) + 0(a 2 s ) (13) 
7r \ m A J 

The one- loop result for the cusp anomalous dimension can be taken from the pioneering paper by Polyakov [6] , 
where it was calculated in the Euclidean space. In Minkowski case, one should change i) — > id to get 

Tlzlr^ocs) = ^C F (tfcothtf- 1). 
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To check the equivalence relation (13), one should just use that Q 2 = 2m 2 (cosh — 1). Taking the 1-loop 
expression for J w (a s ) from [5], 

7™K) = — C F ( r l— \og(w + y/w^T) - l) (14) 

TT V V W 2 - 1 / 

one can also verify, at the one- loop level, the relation (11). 

The equivalence relations (11) and (13) manifest a close connection between the HQET and the theory of 
Wilson lines. The cusp anomalous dimension appears here as a new universal function controlling the IR- induced 
properties of perturbative QCD. All the specifics of the process under study is contained in the dependence of 
r cM sp($,as) on the cusp angle $ determined entirely by the kinematics of the process. 



4. Non-abelian exponentiation theorem 

Using effective fields, the one-dimensional fermions or the effective heavy-quark fields, one can essentially simplify 
a general study of the renormalization properties of the Wilson lines. However, working within the effective field 
theory approach, one can easily miss some properties of the relevant anomalous dimensions which are evident 
if one treats them in the context of the Wilson line formalism. 

For example, in QED, it is possible to get an exact expression for T cusp (d,a). The reason is that, in 
an abelian gauge theory without massless fermions (electrons in QED are treated as massive particles!), the 
vacuum average of a Wilson exponential is an exponential of the vacuum average corresponding to the photon 
propagator: 

(0\Texp(igJ dz'Ariz)) |0) = exp (^f- dz» J^(0\TA^ Zl )A v (z 2 )\0) 
As a result, the cusp anomalous dimension in QED is completely determined by the first loop: 



rg?(«,.)-f*(£) 



There are no higher order corrections. From the IR side, this property is well-known: all the IR singularities 
are given in QED by the exponential of the one-loop result [13]. 

In QCD, the situation is more complicated. First, there arc light quarks u,d and s which can be treated 
as massless at a typical hadronic scale ~ 1 GeV. Second, the gluons are described by non-abelian gauge 
fields. Nevertheless, it is possible to prove a non-abelian exponentiation theorem [17]: the vacuum average of 
a Wilson line operator in QCD is an exponential of some expression, to which, of course, not only the first 
loop contributes. Still, the property observed in QED is partially preserved: the higher-order corrections to 
r C M S p(^,a s ) do not repeat information contained in the one- loop result or, in general, in all the previous loops: 
calculating T cusp {T),a s ), it is sufficient to consider only the diagrams which do not contain the lower-order 
ones as subgraphs, i.e., the essential diagrams are two-particle irreducible in the current channel. As a result, 
the higher-order contributions, in the gluonic sector, are essentially non-abelian. In particular, the second- 
loop term should not contain terms proportional to C F , the square of the one- loop color weight Cp- All such 
terms, present on diagram-by-diagram level, should eventually cancel, and only the terms proportional to the 
essentially non-abelian color weight C f Ca might remain. This observation can be used to check the results of 
the higher-order calculations for the relevant anomalous dimensions. 

The two-loop contribution to T cusp in pure gluodynamics (QCD without massless quarks) was calculated in 
our paper [18]. 2 The massless quark contribution can be also added [9] and the total result is 

r 2 -g»(0,a s ) = (^yC F [-N f ^coth#-l) + C A (±+(^-^y#cothi>-l) 

- cothtf / #Vcoth-0 + coth 2 ?? / dipip{d-ip)cothip (15) 
Jo Jo 



1 . , „„ f , , V'cothV'- 1 , sinh?? 
-sinh2tf / dtp 5 5 — log-r-: — r 

2 J sinh 2 - sinh 2 ip smhip 



independently, a similar calculation was performed by Knauss and Scharnhorst [15]. They failed, however, 
to get rid of all double integrations in their result. 
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It explicitly satisfies the requirement of the non-abelian exponentiation theorem: there are no Cj.-terms in it. 

In the framework of the HQET, the two-loop calculation of the velocity-dependent anomalous dimension was 
performed in a recent paper by Ji [19]. The author presents his result diagram by diagram, in terms of 7 integrals 
ii, . . . ,1?, two of which (ii and I 7 ) are calculated explicitly and 5 others are claimed to be not calculable in 
terms of the elementary functions. The integral Iq even contains a double integration. 3 From our formula 
(15), it is clear that there should be really only three independent integrals containing no double integrations. 
The most worrying is the fact that the Cf-term in [19] is not explicitly equal to zero: it is proportional to 
(/1/3 + I4 — Checking the results of [19], we managed to show that I3 = 0. This observation allowed us 
to simplify another integral, I4. Then, performing a straightforward integration by parts, we obtained that 
^5 = I4. Thus, the total C 2 F coefficient in [19] vanishes, as it should. 

The same statement is valid for the end-point anomalous dimension "f en d(o: s ) 0T i equivalently, for the anoma- 
lous dimension of the effective heavy field in the HQET. In particular, no Cjn-terms really appear at the two-loop 
level. The authors of [16], apparently, have not realized that their two-loop results have this property. 



5. Cusp anomalous dimension for small and large angles 

Incorporating the exponentiation theorem, one can easily obtain an important result mentioned in the introduc- 
tion that the anomalous dimension r cMsp (^, a s ), for large angles i?, is linear in ■& (or log(Q 2 /m 2 )) to all orders 
of perturbation theory: 

T cuap (d,a s )\^ oo =0K(a s ) + O(0°). (16) 

This fact has a simple explanation. In this limit, log(Q 2 /m 2 ) is a typical collincar logarithm resulting from inte- 
gration over the regions where the gluons are almost collinear to v\ or «2, i.e., when they are emitted/absorbed 
at small angle. However, in the two-particle irreducible diagrams (only these really contribute to r cusp (^,a s )) 
there is just one angular integration producing a singularity in the m — > limit. Hence, just a single logarithm 
emerges, and we get eqs.(l), (16). 

Using the explicit expression (15), one can check that our T cusp (i} 7 a s ) is linear in ■& for large However, 
this linearity requirement is not satisfied by the result [19] given by Ji, since it contains an exponentially rising 
contribution 0(exp(i9)) in the term proportional to CfNj. In general, we disagree with his results for the terms 
having this color factor. 

From the expression (2), it follows that the first two coefficients of the perturbative expansion for the 
asymptotic slope K(a s ) are positive. There are good reasons to expect that K(a s ) should be positive to 
all orders. In particular, it can be shown [9] that K(a s ) determines the asymptotic behavior of the parton 
distribution function f(x,Q 2 ) as x — > 1 for large Q 2 . More specifically, the exponent a(Q 2 ) in the asymptotic 
representation f(x,Q 2 ) oc (1 — x) a ^ ) obeys the evolution equation [9] 

<l a(Q 2 ) = K(a s (Q 2 )). 



dQ 

If K(a s ) is a positive function, then a(Q 2 ) is an increasing function of Q 2 , and, for x <~ 1, the parton distributions 
tend to zero faster and faster with increasing Q 2 . And that is just what one is expecting within the QCD parton 
picture: if one probes the structure of a hadron at smaller and smaller distances, the probability to find a parton 
carrying almost all momentum of the hadron should become smaller and smaller. 
In the opposite limit of small cusp angles, eq.(15) gives 



r c „ ap (i?,a s )|0_>o = "d 2 



3tt * \ ir J * \ A \ 54 18 y ' 54 



+ 0{^) (17) 



Apparently, the expansion of the two- loop cusp anomalous dimension (15) contains only even powers of the angle 
■d. In view of the equivalence relation (11) this seems natural: T cusp (d, a s ) depends on $ through w = cosh$, 
which is an even function of However, this simple reasoning is true only if r cusp ($, a s ) is analytic at 1? = 
or, equivalently, if r ) w {a s ) is an analytic function of w in the vicinity of the point w — 1. This fact is not 
immediately obvious, since the explicit 1-loop expression (14) for ^y w (a s ) contains the square root Vw 2 — 1 
capable of generating a cut starting just at w = 1. 



3 We failed to check this integral, since the explicit expression presented in [19] contains a variable u (or a 
function u) not defined in the paper. 



6 



To study the analyticity properties of J w (a s ), it is useful to use its relation to the anomalous dimension of 
the IR factor Fjr (12) and the fact that Fir enters into the factorized expression for the elastic quark form 
factor F(Q 2 , m 2 ). Hence, j w (a s ) possesses only those singularities which the quark form factor has as a function 
of Q 2 = 2m 2 (w — 1). It is well known, that the nearest singularity of the form factor is due to the pair creation 
in the annihilation channel and it corresponds to the point Q 2 = —Am 2 . This implies that ^ w {a s ) is an analytic 
function at w = 1, with the nearest singularity at w = —1. Indeed, using the explicit expressions (3) and (15), 
one can easily check that the bremsstrahlung function has a cut starting at Q 2 = —Am 2 while the two-loop cusp 
anomalous dimension has the only pole singularity at cosln? = — l. 4 

6. Conclusions 

Thus, we demonstrated that three important functions, that appear in different branches of the quantum theory 
of gauge fields, coincide. The same function, the cusp anomalous dimension of a Wilson line, appears when 
one studies the IR singularities of the on-shell form factors or the velocity-dependent anomalous dimension 
in the effective heavy quark field theory. This is because in both cases, all the IR induced terms can be 
factorized from the original amplitudes into the universal IR factor given by the vacuum average of a Wilson 
line. This object, though it is a nonlocal functional of the gauge fields, has rather simple (and well studied) 
renormalization properties. As a consequence, the anomalous dimensions of the Wilson lines govern the IR 
logarithms in the same manner as the anomalous dimensions of the local composite twist-2 operators control 
the collinear singularities one encounters studying the structure functions of deep inelastic scattering. 

In practical aspect, the well-developed machinery of the Wilson loop formalism can provide a good support 
for the heavy quark effective field theory. In particular, we established the equivalence of some recent HQET 
calculations [5, 16] with the earlier results for the Wilson lines [6, 7, 15]. We also corrected the two-loop 
calculation [19] of the velocity-dependent anomalous dimension. 

Acknowledgment. One of us (A.V.R.) is grateful to E. De Rafael and J. Soffcr for the warm hospitality at 
CPT (Marseille). G.P.K. thanks G.Marchesini for useful discussions and support. 
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